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Stable cosmological models driven by a free quantum scalar field 
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In the mathematically rigorous analysis of semiclassical Einstein's equations, the renormalisation 
of the stress-energy tensor plays a crucial role. We address such a topic in the case of a scalar field 
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without an a priori introduced cosmological constant. As a further novel result we shall show that 
these solutions turn out to be stable. 

PACS numbers: 04.62.-|-v, 98.80.Qc, 98.80.Jk 



tained as well a smooth exit from an inflationary phase. 
Since this is a topic partly far away from our goals, we 
shall consider anew such case, namely we study the semi- 
classical Einstein's equation 

Gab ~ 8TTG{Tab)Lu, 

where the left hand side is the standard Einstein's tensor 
whereas the right hand side is the expectation value for 
the stress-energy tensor in the state w. It is a well known 
problem that the latter gives origin to divergences. Hence 
it is compulsory to invoke a renormalisation procedure 
and, within this perspective, we would like to carry on 
our analysis along the lines discussed by Wald, using the 
point splitting regularisation. 

In a series of papers [8, 9], Wald sets out five axioms 
that need to be satisfied in order to have a renormalised 
stress-energy tensor that can be used in order to have 
possible meaningful semiclassical solutions of the Ein- 
stein's equation. Sticking to such a perspective we shall 
show that, in some physically motivated limits, we can 
find a stable solution to the semiclassical Einstein's equa- 
tion. This leads to a great difference from the original 
Starobinsky model where, on the opposite, an unstable 
behaviour is displayed. To this end, we must bear in mind 
the following message already conveyed to us in [T^, |Tl[ : 
the renormalisation of the stress-energy tensor suffers of 
some ambiguities encoded in a modification of the action 
by the addition of terms depending only on the curva- 
ture and on the parameters describing the fields such 
as for example the mass. This arbitrariness is then en- 
coded in the renormalisation parameters present in front 
of this arbitrary terms. In the forthcoming discussion 
we shall fix the renormalisation parameters requiring a 
physically meaningful theory and invoking the principle 
of general local covariance ■ It will also turn out that 
the original result due to Starobinsky in the case of con- 
formal invariant fields corresponds to another choice of 
the renormalisation constants; hence, employing a differ- 
ent criterion, the system under analysis displays a rather 
physically different behaviour. 

For a more mathematically oriented reader a few more 
comments are in due course. Since we are interested in 



I. INTRODUCTION 

A landmark in present days observational cosmology 
has been set by means of the measurement of the type 
lA supernovae red shift which, as a byproduct, proved 
that the Universe is undergoing a phase of accelerated 
expansion. Such a result, also combined with the most 
recent data collected in several other experiments, sug- 
gests that, in order to explain the present state of our 
Universe, we must take into account the presence of a 
"dark energy" playing the role of an effective cosmolog- 
ical constant. From a theoretical point of view we still 
lack a full-fledged satisfactory model for dark energy and 
such a problem was tackled in the past in several ways, 
the most notables being by means either of a yet unob- 
served classical scalar field coupled to wavity P, or of 
a modified theory of gravity itself (see Q and references 
therein for a recent review). 

In the present paper our aim is to consider the backre- 
action of a massive quantum scalar field coupled to grav- 
ity in order to discuss the role played by quantum effects 
in the framework of cosmological models. The interest 
in backreaction effects of quantum fields in cosmology is 
not new since, already in the eighties, Starobinsky Q 
addressed the same topic taking into account a massless 
scalar field conformally coupled to gravity (see also (H). 
The endpoint of Starobinsky seminal paper was the con- 
struction of a graceful exist from a de Sitter phase of 
rapid expansion. Using quantum property of the source 
fields he observed that such a de Sitter spacetime is an 
unstable solution of the semiclassical Einstein's equations 
(see also [6] ) . More recently, in ^3] , Shapiro and Sola also 
considered the massive case in a similar way. They ob- 
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solutions of the semiclassical Einstein's equation, where 
quantum matter acts as a source for the gravitational 
field, we need to employ a quantisation scheme indepen- 
dent from the spacetime itself. Such a conceptual prob- 
lem was recently addressed in a work due to Brunetti 
Fredenhagen and Verch • They showed that it is pos- 
sible to simultaneously quantise on all spacetimes and the 
quantisation scheme in this framework corresponds to as- 
sign a functor between the category of spacetimes (Man) 
and the category of local Algebras Loc generated fields. 
Furthermore such a functor transforms covariantly under 
any local transformation. Unfortunately, while also fields 
transform covariantly under isometrics, a similar conclu- 
sion cannot be drawn for states. Therefore, since we are 
interested in expectation values of fields, we are forced to 
select a class of the mentioned states enjoying some suit- 
able physical properties and in the framework of FRW 
spacetimes, this naturally leads to select the class of the 
so-called adiabatic states. Starting from these premises 
we are now ready to use, within this abstract scheme 
of analysis, quantum matter as a source for the grav- 
ity whereas the role of Einstein's equations will select 
a particular set of objects in Man, as a sort of consis- 
tency check. To rephrase, even if we can quantise in all 
the spacetimes simultaneously, once a family of states is 
chosen, only in few of those spacetimes the semiclassical 
Einstein's equations hold true. 

After fixing some notation, in the next section we shall 
recall briefly the renormalisation procedure we shall em- 
ploy. In the third section we shall perform a suitable 
choice for the quantum state and, then we will discuss 
the associated solutions of the semiclassical Einstein's 
equations. In the fourth section we shall justify this hy- 
potheses by means of physical motivations. Finally some 
conclusion are drawn in the last section. 



A. Einstein's equation and cosmological 
backgrounds 

To set notations and conventions, let us clarify that 
our aim is to consider spacetimes whose metric is used in 
the description of the Universe. Hence, we stick to the 
standard convention of requiring the Cosmological Prin- 
ciple to hold true; this straightforwardly leads to the full 
class of Friedmann-Robertson- Walker metrics and, par- 
ticularly, here we shall only consider those with a spatial 
flat section. In a Cartesian reference frame the metric 
reads 

ds^ = -dt^ + aitfSijdx'dx^, i,j = 1, ...,3 (1) 

where a(t) can be interpreted as usual as the expansion 
factor and it is the only function to be determined out 
of (the semiclassical) Einstein's equations. A standard 
calculation shows that these can be reduced to an identity 
at a level of traces i.e. 

-R = 8n{T)^, (2) 



together with the conservation law for the stress-energy 
tensor, namely 

v"(r,b)^ = o. (3) 

As already remarked in the introduction, {T)^^ stands 
for the expectation value of the stress-energy tensor. We 
shall deal with this issue more in detail in the forthcoming 
discussion. As a last remark we wish to recall that ^ 
and ^ are actually not equivalent to a single but only 
to a set of Einstein's equations which differ only by a 
conserved traceless tensor T^^; such arbitrariness is fixed 
upon imposing suitable initial conditions. 

II. MASSIVE SCALAR FIELD. 

As we already emphasised in the introduction, we shall 
employ a real scalar field (j) as the prototype to discuss 
the quantum behaviour of classical matter on a FRW 
background ([!]). Therefore the classical dynamic of our 
system is governed by 

P^^O, P -.^ -D + ^R + m^, (4) 

where ^ G M, i? is the scalar curvature whereas m is 
the mass of the field. Bearing in mind that, unless 
stated otherwise, our convention for the metric signa- 
ture is (— , -f , -f , -f ), (P) entails the following identity 

R — 6 ( - + % ) where each dot stands for derivation 

with respect to t. In the next we shall indicate H = d/a. 
Setting S, = ^ corresponds to the so-called conformal cou- 
pling. 

A. Quantisation procedure: States and Hadamard 
condition 

In this paragraph we shall start dealing with the quan- 
tum behaviour of the solutions of fl]) and, to this avail, 
we shall stick to the realm of the algebraic formulation 
of quantum field theory. Since a detailed analysis of the 
main ingredients and results would require a review on 
its own just for the massive scalar field, we shall point 
an interested reader to [H, Therefore, to cut a 

long story short, let us state that, to our purposes, it 
suffices to remember, that being the FRW spacetime, 
globally hyperbolic, it exists a standard procedure to 
assign a *— algebra, say W, out of ^ ^3 After- 
wards we need to add a further ingredient, namely a state 
uj : W ^ C, which is the key tool out of which we can 
calculate the relevant objects i.e. expectation values of 
the fields on that state, more commonly referred to as 
n— point functions which we shall denote from now on as 
= {4'{^i) ■ ■ ■ 4'{xn))- From a formal perspective these 
objects must be thought as distributions in 'D'{M'^) and 
the singular structure, proper in general of distributions, 
arises whenever we perform in aj„ a coincidence limit. 
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Therefore, in order either to formulate a mathemati- 
cahy meaningful field theory either to construct a theory 
which allows us to perform calculations going beyond the 
pure formal level, the selection of a suitable class of states 
is one of the main, if not the most important, task. To 
this avail, we shall impose some reasonable constraints 
and the first requires us to restrict the attention to the 
so-called quasi-free states. These are characterised by the 
following property: all the odd n— points functions van- 
ish while all the even can be reconstructed out of sums 
of products of the two-points function. In other words 
quasi free states are fully determined once 0^2(2;,?/) is 
known. In the forthcoming sections we shall display how 
the above requirement is relevant to our discussion. In 
particular we shall show that also the stress-energy ten- 
sor can be fully determined only out of W2 and this is 
the key non-geometrical ingredient in the semi-classical 
Einstein's equation. 

Nonetheless "quasi-free" is not a sufficient requirement 
for our Lo to satisfy and, particularly, a second and most 
important hypothesis must be imposed, namely the state 
shall be Hadamard. On a practical ground, from such a 
condition we can infer that the singular structure for the 
two-points function is fixed as 

ax,y 

where a is half of the square of the geodesic distance in 
the FRW background. The functions u, v and w, also 
known as Hadamard coefficients, are smooth and u, v 
can be uniquely determined once the equation of motion 
and the metric of the underlying background are fixed. 
In the above expression it turns out u is the square root 
of the so-called van Vleeck-Morette determinant which 
depends only on gab, i-e. u can be reconstructed only 
out of the geometric properties of the manifold on which 
our fields live. On the opposite, w is the contribution 
to the Hadamard function which depends upon the state 
we have selected. Therefore all the information of the 
singular part in (O is encoded in 

TT, \ 1 / uix, y) / X , / x\ 

Hy^^ v)^ —( — T + "(^' y) log ^0^^' v) ' 

Stt^ \(T[x,y) J 

which has a universal structure in every Hadamard state. 
Hence this is the contribution that we can subtract from 
the two-points function in order to get a smooth be- 
haviour; in other words this amounts to regularise the 
state. As a notational convention, from now on, we shall 
refer to v{x, x) by means of the symbol [v\. Furthermore 
v{x,y) admits an asymptotic expansion in powers of the 

geodesic distance: v{x,y) — ^ Vn{x,y)a"{x,y). In the 

n=0 

forthcoming discussion the coefficient vi will play a dis- 
tinguished role. 



B. Stress-energy tensor 

The stress-energy tensor for a quantum real scalar field 
(/> with mass m and coupling to curvature ^ can be written 
as 

Tab := da4>db^ - ^9ab {dc^d^^ + m'cf') - ^Vadb4>^+ 

H {Rab - l^ab) 0' + - ^ gabOcj)^. 

Since the key ingredient to our analysis is the trace and 
the conservation equation for Tab, let us switch from the 
previous formula to 

T = -3 Q - ' V,T% = 0. 

We stress to the reader that, here, we employ a non- 
standard form for Tab, i-e. it differs from the more famil- 
iar one by a term proportional to ^ {{P(j))(j) + (j){P(j))) gab 
[l5l |. At a classical level this contribution vanishes since, 
on shell, = 0, but nonetheless it represents an im- 
portant feature in a fuU-fiedged analysis of the under- 
lying quantum theory, since, in this case, it is different 
from zero. Furthermore, encompassing such a term in 
the stress-energy tensor, automatically accounts for the 
trace anomaly which, on the opposite, was usually added 
by hand. As shown in [l^, [HI, [H, this automati- 
cally arises in the quantum theory once the point splitting 
regularisation is performed. We also exploit the latter 
in order to regularise the operator Tab in order, subse- 
quently, to calculate its expectation value on a quasi-free 
Hadamard state. Such an expression would be quite cum- 
bersome in the text and also of little avail; therefore an 
interested reader can refer to the appendix A.l for more 
details. 

Notice that the envisaged conservation equation for 
the quantum stress-energy tensor, namely ya{T°'^)ui = 0, 
holds true due to the following identities 

87r2(0P0)^ = 6[t;i], Sn"" {{V a(t>){P4>))^ - 2Va[wi], 

where [vi] is here explicitly given in the appendix in for- 
mula (|A1|) . The heritage of such a conservation law is 
the change of the expectation value for the trace of Tab 
by means of a purely quantum term: 

<r,..(-3(l-,)n-„.)H.M, 

where the dependence upon the state is encoded in the 
term [w]. 

To conclude, we point out to a potential reader that, 
due to [vi], the above trace is non vanishing also in a 
conformal field theory 

C. Remaining freedom in the definition of Tab 

By means of point splitting regularisation we have fixed 
the expectation value of {T)i^ in the so-called minimal 
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regularisation prescription, namely we have only sub- 
tracted the singular part form the two-points function. 
Nonetheless, as discussed by Wald 0, in the renormal- 
isation prescription, there is still a freedom of geomet- 
ric nature. In detail we can add a tensor tab written 
only in term of the local metric and such that it satisfies 
V^iah = without either affecting the equations of mo- 
tion for the matter either violating the first four axioms 
introduced and discussed in Wald paper. The conser- 
vation equation for tab is not the unique constraint we 
may wish to impose on such a tensor and, in particular, 
a further natural requirement would be that tab behaves 
as Tab under scale transformations. In other words this 
implies that tab arises out of the following variation 



tab — 



ah 



A^R" + B^RabR''\ 



being A and B just arbitrary real numbers. Leaving 
the details of the above construction and analysis to 
0, [13, [iH , we shall only stress that the trace of tab turns 
out to be proportional to Di? independently from the 
choice of A and B. This is an unavoidable arbitrariness 
in the employed scheme and, as a byproduct, it leads us 
to think of A and B as renormalisation constants on their 
own. We are now able to compute the trace of the whole 
quantum modified stress-energy tensor: 



87r2 



2H 
87r2 



-cDi?, 



where c is a linear combination of A and B and it rep- 
resents the freedom in the renormalisation procedure we 
exploited. Eventually, c will be chosen in order for the 
trace to satisfy the requirement coming out of the fifth 
Wald's axiom (still see Q); in other words there must 
be no derivatives of the metric with degree higher than 
2 in the expectation values of Tab- The remaining renor- 
malisation ambiguity is encoded in the expectation value 
of the filed {(f>'^)uj; we shall come back later to this point 
fixing the ambiguity by physical motivation. 

We stress that a similar observation brought interest 
in the so-called modified theory of gravity also known 
as f{R) gravity. Nonetheless the view we wish to push 
home is the following: adding tab does not come from a 
modified gravitational action, but it only originates form 
the employed renormalisation scheme, i.e. it must be an 
effect coming from quantum matter. Naturally this does 
not exclude that such a perspective cannot provide hints 
on how a candidate theory of quantum gravity interacts 
with quantum matter. As a final comment we would like 
to stress that the above is the subtlest point in the whole 
construction. We used an expression for the stress-energy 
tensor which is suitable in order to deal with semiclassical 
Einstein's equation. Nonetheless such a modification is 
not artificial, corresponding as a matter of facts just to 
a specific choice of the renormalisation constants arising 
out of the employed scheme. 



III. EVOLUTION EQUATION OF THE MODEL 

In the case of conformal coupling ^ — 1/6, equation 
written in terms oi H ^ a/a, becomes 



6{H + 2H^ = -SnGm 



30 



G 



HH^ + H 



m 
T 



(6) 



The aim of this section is to analyse in detail the pos- 
sible solutions of (O under some specific hypotheses on 
the expectation value for {4>^)ui- Particularly we shall 
show that a de Sitter space with a specific curvature will 
appear as a stable solution. 



A. Conformal invariant case: stability of de Sitter 
phase 

As a starting point we shall deal with the scenario in 
which TO = 0, already encompassed in Starobinsky paper 
[1] (see also [5]). As remarked above, there is no need to 
select a specific state and an ordinary differential equa- 
tion rules the evolution of H. Hence, setting to = in 
([6|), we end up with 



(7) 



Here Hq = i^^: depends on the Newton constant and it 
has an order of magnitude of 24 times the inverse Planck 
time. Let us notice that, out of the right hand side of (O, 
we can extract two critical points; therefore ([7]) admits 
two constant solutions, namely H{t) — and H{t) = 
— corresponding respectively to a Minkowski 

spacetime and to a de Sitter one. Suppose now to assign 
an initial condition at a fixed time to such that H{to) ^ 
and H{to) ^ 77+; we are interested to realize if the 
solution interpolating such an initial condition flows at 
large times either to or to Hj^ i.e., in order words, 
whether these two critical points are stable or not. To 
bring such task to a good end, we simply need to notice 
that ([7]) is integrable as: 



Ke 



4t _ g2/_ff 



H + H. 



H - H4 



1/H4 



(8) 



where K stands for the integration constant to be fixed 
out of the initial condition Hq. Depending on such last 
value, all the solutions H{t) flow either to or to 
Hence both critical points turn out to be stable. This re- 
sult is different from the classical outcome of the analysis 
due to Starobinsky [4] (see also Vilenkin and Ford [5, 6]). 
The price to pay, in order to achieve such a result, is a 
choice by hand of a renormalisation constant. It turns 
out to be an addition of a tensor written only in terms of 
the metric and such operation introduces in the theory 
a scale- length, as already discussed by Wald in 0. We 
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have to stress that, on the dark side, the above de Sit- 
ter solution cannot describe the present days form of the 
universe being ~ 6.4 x lO^'^s^^ i.e. many orders of 
magnitude bigger then the present measured Hubble con- 
stant (2.6 ±0.2) X lO^i^s"^ On the bright side, instead, 
we have shown that, encompassing the full quantum ef- 
fects, we are lead to find a stable de Sitter solution even 
if no cosmological constant is present in the equations. 



B. Massive case with ^ = 1/6: stability of the de 
Sitter phase, effective cosmological constant 

In this section we switch from the massless to the mas- 
sive case. The most important difference is the follow- 
ing: the righthand side of ^ depends explicitly upon 
the state via the expectation value of cjp. The expec- 
tation value of {4>^)ui on a general Hadamard state uj is 
-|- arn? -\- PR, where a and f3 are renormalisation con- 
stants encoding the ambiguities still present in the proce- 
dure. We assume for the moment the existence of a set of 
Hadamard states one for each spacetime whose metric 
is of the form ^ being H — a/a and {(jy^)^ — am? + (3R. 
We shall see later that this assumption turns out to be 
an approximation of the expectation values of the fields 
computed on the adiabatic states of FRW in the limit 
where m? >> R and m» H. Moreover, by the princi- 
ple of general local covariance [HI [13 , we are entitled 
to fix the renormalisation constants once and in the same 
way for every spacetime we are considering. Then the ex- 
pectation value of (0^)2 on the states we are considering 
takes the following values: 



= arr? + PR, 



(9) 



on all the considered FRW spacetimes. Therefore, taking 
into account these remarks, ([6]) takes the following form: 

H - H^) = -H^ + 2HIH^ + M, (10) 

where and M are the following two constants with 
the following values 



2 180^ 



STT^lSOm^/?, 
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As in the previous section, the right hand side of (fTO|) 
displays at most two critical points amounting to 



Hi = 2 ± 



M, 



(11) 



both corresponding either to a de Sitter phase or to a 
Minkowski phase. A straightforward analysis shows that 
both H{t) = H± appear to be stable since all the so- 
lutions flow to either one of the two fixed points. It is 
remarkable that the existence and the stability behaviour 
of the latter is left unchanged whether the right hand side 
of ^ is modified adding a term such as Aa~^{t), being 
A e M and A a constant of suitable dimension. It is 



also interesting to notice that a formula similar to (fTTj) 
already appeared in although, in the cited paper, a 
classical cosmological constant has been introduced from 
the beginning. At this stage our simple model depends 
on three parameters a, P, m. A minimal and, to a certain 
extent, compulsory choice is to require Minkowski as a 
solution of our system. This amounts to fix a = (327r^)~^ 
which, on the other hand, entails M = 0. The form of 
the solution is then equal to that of the massless case ([8]) , 
where one of the fixed points corresponds to a Minkowski 
space - H{t) = -, while the other fixed point H{t) — 
corresponds to de Sitter. With respect to the massless 
conformal factor, here we can fine-tune the parameters P 
and m in such a way for iJ+ to be small enough in order 
to account for the present measured value of the Hubble 
constant. Hence, heuristically speaking, our system be- 
haves as if an effective cosmological constant enters the 
fray without even being present at the beginning and this 
is a strict consequence of encompassing the full quantum 
properties of the field. As a further remark we would 
like to notice that ^ displays, for a large class of initial 
conditions, an early time phase of rapid expansion which 
is a prerequisite feature of modern models for studying 
the early stages of evolution of the Universe. This is 
in sharp contrast with the canonical paradigm according 
to which quantum effects should account only for small 
fluctuations with respect to the classical behaviour. On 
the opposite, even in the most simple example of a mas- 
sive scalar field and with the most simple assumptions, 
our system displays a behaviour which drastically differs 
from the one we could a priori expect only from a clas- 
sical analysis. Hence this suggests that, when dealing 
with scalar fields on a FRW background, one should al- 
ways perform a full-fledged analysis of the semiclassical 
behaviour of the system since the quantum contributions 
appear to be hardly negligible as one can also infer from 
flgure [TJ 

As a flnal comment we would like to stress that, in a 
neighbourhood oi H = the found solution ([5]) looks 
rather similar to the one of a classical flat universe with 
cosmological constant filled with radiation. As a matter 
of fact, in that case H{t) = Atanh(2(t — ^o)^) where A 
is a constant related to the cosmological constant, and it 
can be inverted as 



Ke 



it 



H + A 



H-A 



l/A 



which looks very similar to ^ when H ^ A and = A; 
this corresponds to the dashed line in figure[T] The quan- 
tum effects are not important only around H = where 
([8]) looks like H{t) in a flat universe filled only with radi- 
ation, namely the dotted line in the flgure [TJ Eventually 
we would like to stress that considering the upper brunch 
of the solution, in the past, it displays the behaviour of 
a classical flat universe with a kind of matter such that 
p = A a(t)'^'^ . Even in this regime quantum effects are 
not negligible. As a further remarkable consequence of 
the analytic form of H{t), it turns out that the singular- 
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FIG. 1: Here the dashed hne corresponds to the behaviour of 
as a function of time t (normalised with respect to l/_ff+) 
in a FRW universe with a non vanishing cosmological constant 
and filled with radiation, while the dotted line stands for the 
lone classical contribution coming from radiation. Conversely 
the continuous line depicts the outcome of our model where 
quantum effects are also taken into account. 



ity at i = to coincides with null past infinity in the flat 
spacetime conformally related to ((1]); hence it descends 
that the particle horizon is not present. Therefore any 
pair of points in the underlying background was casually 
related in the past, and, thus, as a byproduct, such prop- 
erty of our model could provide a solution to the problem 
of homogeneity. 



IV. 



EXPECTATION VALUE OF 4? ON THE 
ADIABATIC VACUUM 



In the preceding section we have seen that, assuming 
a suitable form of (c^?)., two stable de Sitter phases can 
arise as solutions of the semiclassical Einstein's equation. 
We would like to give a justification for our assumption, 
namely we shall show that there is a regime in which it 
is is valid. Here we restrict our attention to the case of 
a massive scalar field with a conformal coupling to the 
metric. The first observation is that, if we select the 
Bunch-Davies state ujb [3 ^ Sitter spacetime and 
if we compute the renormalised version of the expectation 
value of 0^, we obtain a constant that depends only on 
the mass m and on H . With this observation we can 
immediately conclude that the two fixed points Hit) — 
and H{t) = discussed above are really exact 
solutions of the semiclassical Einstein's equation. In the 
next we shall select a class of states that, in the limit of 
a large mass, shows an expectation value for {(jP)uj that 
is of the type am^ + (3R. 



A. Adiabatic states and large mass expansion 

We would like to select here the class of adiabatic 
states, i.e. those introduced by Parker fl^ in order to 



minimise particle creation (see also [20| for a derivation 
of the expectation values of the stress tensor). Much 
work has been done also recently in order to make the 
definition of these states precise [2ll [2^ . [2^. In order 
to write the two-points function of these states we follow 
the construction as in Parker [l9[. In the case of confor- 
mal coupling it is convenient to use the conformal time 

Therefore the two-points 



T defined as r — rn = C -^ttt 

" Jto a{t) 

function of such kind of states is 



LU{X1,X2) 



1 



1 



87r3 a{Ti)a{T2) 



ik-(xi-X2) 



above four vectors and X; are three vectors 

whereas |k| stands for the length of the spatial vector 
k. The functions ^'fc(r) are solutions of a differential 
equation with a suitable normalisation condition: 



'^kiT) — ^k{T) =i. 



Each \E'fc(T) can alternatively be written in the following 
way: 



1 



In the adiabatic approximation r2fe(T) is a function con- 
structed recursively in the following way: 



and 



(n+l)^ 



(t) = k'+m'a{TY+- 



3(nr{r)\ If^i") (r) 



2 ^':Hr) ' 
(12) 

where the prime stands for the derivation with respect 
to T. The n-th order approximation consists then in the 
substitution of flk with fi^"^ in ^'^(t) and we shall indi- 

(n) 

cate with the counterpart for the two-points func- 
tion of the state. Nonetheless one should bear in mind 
that this recursive procedure does not have nice conver- 
gence properties though, thanks to the work of Junker 
and Schrohe [13], we know that the state constructed 
in this way is an adiabatic state in the sense that Wj"' 
have a certain Sobolev wavefront set. Hence, if n is large 
enough, we can use the approximated state in order to 
build the stress-energy tensor or the expectation value 
of (jP . In particular, we can compute the approximated 

expectation value {4>^)in) = lini (tj^n"^ (^i 2/) ~ H{x,y)), 
which, more explicitly, becomes 



1 



I in) 



47r2 a(r)2 



dk 



(0), 
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Above a' and fi' need to be interpreted as renormalisation 
constants. An exact computation of this integral can be 
very difficult to perforin, hence we will show only how 
to compute an expectation value in the limit of a large 
mass, namely, assuming that H{t) is a smooth function 
and >> R. In this case, if furthermore n > 2, it is 
possible to expand the integral in powers of l/m?, as: 



/(«) 



am^ +PR + 



where a and P are slightly different from the one written 
before. In the large mass limit we shall simply consider 
{4''^){n) — ct^^ + PR- The result should be read as a 
confirmation for the approximation we have done in the 
preceding section. 



INTERPRETATION OF THE RESULTS AND 
FINAL COMMENTS 



Shapiro and A. A. Starobinsky for useful comments and 
remarks. 

APPENDIX A: POINT SPLITTING 
REGULARISATION OF THE STRESS-ENERGY 
TENSOR 

Let UJ2 be the two-points function of a quasi free 
Hadamard state. The expectation value of the stress- 
energy tensor regularised by means of the point splitting 
procedure is: 



{Tab)uj{z) 



lim 



-2a^adb+dadl)+^[Rab{z) 



R{z) 



lab 
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{uJ2{y,x) - H{y,x) +uj2ix,y) - H{x,y)) 



In the present paper we have shown that, when dealing 
with cosmological models, quantum effects are not negli- 
gible even when we consider basic models. As a matter of 
facts, our analysis displays that, from a careful analysis of 
the expectation values of the renormalised stress-energy 
tensor, it arises an effective cosmological constant which 
can be interpreted as dark energy. 

Such a feature is manifest if we take into account 
a massive scalar field propagating in a curved back- 
ground, although we envisage that similar effects would 
be present if we consider other kinds of fields. Further- 
more we have seen that a de Sitter solution appears as a 
stable fixed point of the semiclassical Einstein's equation 
and, to a certain extent, also a phase of rapid expansion 
can be foreseen in the model. We also believe that, since 
the found results, and particularly the stability of the 
de Sitter solution, are based upon a modification of the 
point splitting procedure by a pure gravitational term, 
this could be read as an hint for future study of quan- 
tum gravitational models interacting with matter. To 
this avail it also seems interesting to pinpoint that, even 
considering the one-loop corrections to the action of an 
f{R) theory, one is lead to a stable of de Sitter solution 
[2J, [2^. Furthermore, also in this last case, stability is a 
joint effect of quantum theory and classical gravity and 
this is a behaviour which a lone f{R) = R^ term does 
not display. 
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where the prime stands for a derivative in y whereas the 
one without prime is a derivative with respect to x. A 
reader should notice that, in the last part of the equa- 
tion, there is a symmetrisation done at the level of two- 
points function and that H{x, y) is the singular part of 
the Hadamard series. 



1. [vi\ coefficient in the cosmological case 

Since it is a relevant datum in our procedure, we pro- 
vide the explicit expression for 2[vi\ = [a2]/2, being 02 
the^chwinger-de Witt coefficient as derived at pag. 194 
in [16] with the choice of = ^i? + m?, (see also |26| ) 
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Furthermore, assuming that the metric has the form of 
a flat FRW universe ((T|) and writing H — a/a, [vi] takes 
the following form 



2W = -1(hH^ + H.) 



12 V 5 
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